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ABSTRACT 


The  effects  of  a bandpass-limiting  satellite  channel,  such  as  a B=25-kHz 
hard-limiting  FLEETSAT  channel,  on  several  modulation  techniques  were  simulated 
The  modulations  included  BPSK  and  various  forms  of  QPSK,  including  offset 
QP5K  and  continuous-phase  minimum  frequency  shift  keying  (MSK).  Coherent 
matched-filter  receivers  and  both  uplink  and  downlink  noise  were  assumed. 

For  a bit  error  probability  of  0.005  and  a SNR  loss  of  less  than  1 dB 
caused  by  the  nonlinearity  in  the  FLEETSAT  case,  the  maximum  standard  data 
rate  R was  19.2  kbps  for  BPSK  but  32  kbps  for  offset  QPSK,  MSK  and  QPSK. 

Rates  of  24  kbps  for  BPSK  and  48  kbps  for  offset  QPSK  and  MSK  were  attainable 
at  this  same  reliability  with  a SNR  degradation  of  about  4 dB  and  2.5  dB, 
respectively.  Larger  degradations  are  expected  at  lower  error  rates. 
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I.  INTRODUCTION  AND  SUMMARY 

In  this  section  we  review  the  motivation  for  this  study  and  present  the 
principal  results  of  our  simulations.  The  simulation  model  is  discussed  in 
Section  II,  and  the  performance  is  analyzed  in  Section  III. 

1 . 1 Background 

The  next  generation  of  military  communications  satellites  denoted  here  by 

•f 

FLEETSAT’s  will  have  a number  of  25-kHz  wide,  bandpass  hard-limiting  channels 

at  VHF/UHF  frequencies  [1].  When  contemplating  the  most  effective  use  of  these 

* 

resources,  one  of  the  issues  that  arises  is  the  maximum  data  rate  feasible 
through  each  channel  for  a given  modulation  scheme.  This  issue  is  addressed 
here  for  coherent  matched-filter  receivers  and  modulation  schemes  employing 
various  forms  of  antipodal  signaling.  The  major  modulations  are:  binary 
phase  shift  keying  (BPSK) , offset  or  staggered  quadriphase  shift  keying  (SQPSK) , 
and  minimum  shift  keying  (MSK) , a form  of  continuous-phase  binary  frequency  shift 
keying  (CPBFSK) . The  main  objectives  are  to  determine  the  better  modulation 
techniques  and  to  estimate  performance  losses  with  respect  to  a linear  wideband 
channel  as  a function  of  data  rate.  The  important  issue  of  how  the  distortion 
introduced  by  the  narrowband  channel  affects  the  synchronization  techniques  used 
with  each  modulation  scheme  for  acquiring  and  tracking  both  symbol  timing  and 
phase  or  frequency  offsets  will  not  be  treated  here. 

We  model  the  satellite  as  a narrowband  filter  followed  by  an  ideal  bandpass 
limiter  that  removes  the  envelope  modulation  of  its  input  (see  Fig.  1.1).  Since 

ri 

Preferably  selected  from  the  standard  families  75  x 2 bps,  n = 0,  1,  2,..., 

and  8000  x N bps,  N = 1,  2,  3, 

1* 

FLEETSAT’s  will  also  have  5-kHz  channels  and  a 500-kHz  channel;  the  results 
of  this  report  apply  to  these  channels  with  appropriate  bandwidth  normalization. 
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Fig.  1.1.  Computer  simulation  model. 


the  space  segment  of  the  system  is  well  past  the  design  stage,  we  view  the  sat- 
ellite as  fixed  and  have  concentrated  our  effort  in  this  area  on  developing  a 
reasonable  model  of  a satellite  filter  from  data  made  available  to  us  through 
the  Navy. 

Typical  characteristics  of  a FLEETSAT  filter  are  shown  in  Fig.  1.2.  The 
filtering  is  actually  done  by  a cascade  of  two  filters  denoted  by  A and  13  with 
a combined  3-dB  bandwidth  closer  to  27  kHz  than  25  kHz.  The  actual  center  fre- 
quency of  this  filter  may  vary  by  + 1 kHz  according  to  the  FLEETSAT  satellite 
specification.  In  this  report  we  normalized  bandwidths  and  data  rates  to  a 
3-dB  filter  bandwidth  of  25  kHz  to  be  slightly  conservative. 

The  phase  deviation  curve  of  Fig.  1.2b  covers  most  of  the  filter* s 3-dB 
bandwidth.  The  FLEETSAT  specif cation  called  for  near  linearity  over  the  center 
80%  of  the  useful  bandwidth.  The  linear-phase  reference  evidently  resulted 
from  a best-fit  straight  line  through  the  band  center.  A slightly  smaller 
slope  tangent  to  the  phase  characteristic  at  the  zero  frequency  deviation  would 
result  in  a phase  deviation  of  about  30  degrees  at  10  kHz  and  -6  degrees  at 
-10  kHz. 

We  have  simulated  several  filters  to  obtain  information  on  system  perfor- 
mance sensitivity  to  a difference  in  spectral  roll-off,  phase  nonlinearity  and 
center  frequency  offset.  This  might  suggest  new  specification  criteria  in  the 
future . 

Although  the  space  segment  is  considered  inaccessible,  we  may  well  have 
the  freedom  to  modify  communication  terminals  belonging  to  the  earth  segment. 
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a.  Response  of  filter  A (filter  _B?s  response  is  similar). 


Fig.  1.2.  Typical  characteristics  of  a FLEETSAT  channel  filter. 
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b.  Phase  deviation  relative  to  linear  reference. 

Fig.  1.2.  Continued. 


5 


This  is  the  principal  reason  for  considering  modulation  techniques  in  addition 
to  the  more  common  BPSK  used  with  the  WSC-3  modem  [2],  for  example.  A next 
step  would  be  to  investigate  sophisticated  (possibly  adaptive)  receiver  struc- 
tures that  are  better  matched  to  the  distorted  waveforms  transmitted  by  the 
nonlinear  narrowband  satellite. 

1. 2 Conclusions 

Six  modulation  techniques  that  yield  equivalent  coherent,  matched-f liter 
performance  in  the  presence  of  additive  white  Gaussian  noise  (AWGN)  were 
simulated  for  the  narrow  bandpass-limited  satellite  model.  Of  these  modula- 
tion techniques,  BPSK  performed  far  worse  than  the  others  at  the  same  data 
rates.  MSK  and  SQPSK  proved  to  be  the  best  modulations  overall,  followed  closely 
by  quadriphase  shift  keying  (QPSK) , and  more  distantly  by  a modulation  scheme  we 
call  alternating  QPSK  (AQPSK)  and  a form  of  continuous-phase  quadrif requency 
shift  keying  (CPQFSK)  (see  Section  2.2.).  In  this  section  we  shall  focus  mainly 
on  the  best  modulations  MSK  and  SQPSK  because  they  are  so  attractive,  and  the 
worst  modulation  scheme  BPSK  because  it  is  so  common. 

The  effective  SNR  losses  for  equal  uplink  and  downlink  SNR’s  compared  to  an 
AWGN  channel  at  a bit  error  probability  of  P^  ~ 0.005  are  depicted  in  Fig.  1.3 
for  BPSK  and  MSK.  SQPSK  performance  is  not  significantly  different  from  the  MSK 
results  of  Fig.  1.3b.  The  data  represent  most  of  the  simulations  performed  for 
these  modulations  using  different  satellite  filter  models. 

The  number  beside  each  datum  gives  the  ratio  of  the  filter  3-dB  bandwidth  B 
and  the  data  rate  R for  that  simulation.  The  numbers  in  parenthesis  for 
filter  6 in  Fig.  1.3  refer  to  phase  deviations  defined  below.  The  solid  curves 
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FASTER  SPECTRAL  ROLL-OFF  — ► SLOWER 


a.  BPSK. 

Fig.  1.3.  Performance  loss  vs.  filter  bandwidth  ratio  for  various  3-dB  bandwidth 
to  data  rate  (B/R)  ratios  (equal  uplink  and  downlink  SNR's). 
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b.  MSK . 

Fig.  1.3.  Continued. 
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are  rather  tenuous  estimates  of  performance  loss  versus  filter  selectivity 
for  constant  B/R  ratios.  The  Xfs  mark  the  best  FLEETSAT  loss  estimates.  Al- 
though these  estimates  are  not  very  sensitive  to  small  changes  in  P^,  the 
losses  would  increase  markedly  with  a decade  or  more  decrease  in  P^. 

The  performance  degradations  are  plotted  as  a function  of  the  ratio  of 
the  filter  20-dB  and  3-dB  bandwidths.  In  each  part  of  Fig.  1.3,  six  labeled 
vertical  lines  mark  these  ratios  for  the  filters  simulated,  while  the  seventh 
vertical  line  labeled  FLEETSAT  denotes  the  bandwidth  ratio  for  the  cascaded 
filter  of  Fig.  1.2,  our  current  information  as  to  the  satellite  filter  charac- 
teristics. 

This  measure  of  spectral  roll-off  for  a filter,  the  20-dB  to  3-dB 
bandwidth  ratio,  is  also  used  in  [3].  This  ratio  happens  to  vary  from  -7% 
to  +13%  for  the  filters  simulated,  relative  to  the  FLEETSAT  filter.  The 
purpose  of  presenting  the  data  in  this  form  is  to  show  how  the  performance 
loss  depends  on  filter  selectivity,  defined  as  the  above  bandwidth  ratio  in 
this  instance.  This  enables  us  to  obtain  a more  accurate  estimate  of  the 
true  degradation  for  the  FLEETSAT  filter,  once  differences  in  phase  non- 
linearities  are  taken  into  account. 

Two  versions  of  filter  6 (see  Fig.  1.3)  were  employed,  one  with  a nominal 
phase  deviation  of  12  degrees  at  10  kHz,  and  the  other  with  a deviation  of 
36  degrees.  From  the  elliptic,  Chebyshev  and  Butterworth  phase  characteristics 
in  Appendix  B,  it  can  be  verified  that  the  Butterworth  and  Chebyshev  phase 
deviations  are  roughly  comparable  in  magnitude  to  the  12-degree  and  36-degree 
deviations  of  filter  6,  respectively;  the  elliptic  phase  deviation  is  an 
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order  of  magnitude  larger. 

On  intuitive  grounds  and  from  the^  data  of  Fig.  1.3,  we  claim  that  the 
larger  the  bandwidth  ratio,  the  more  sensitive  the  performance  loss  is  to  a 
change  in  phase  nonlinearity.  This  is  because  the  magnitude  characteristic 
dominates  the  nonlinear  phase  effects  in  regions  of  large  attenuation.  Thus, 
since  the  12-degree  filter  6 more  closely  approximates  the  known  phase  devia- 
tion for  the  FLEETSAT  filter,  more  weight  should  be  given  to  the  12-degree 
data  than  the  36-degree  filter-6  data. 

In  drawing  the  curves  of  Fig.  1.3,  then,  the  data  for  filter  5^  the 
Butterworth  filter  and  the  12-degree  filter  6 were  weighted  more  heavily, 
followed  by  the  Chebyshev  and  36-degree  filter-6  data.  The  filter-3  data 
were  given  less  emphasis  because  the  phase  characteristic  is  linear.  On 
the  other  hand,  the  phase  of  the  elliptic  filter  is  too  nonlinear.  The  inter- 
sections of  the  degradation  curves  with  the  FLEETSAT  lines  of  Figs.  1.3a  and 
b indicate  our  estimates  of  the  performance  losses  to  be  anticipated  for  the 
poorest  modulation  scheme  BPSK  and  the  best  modulations  MSK  and  SQPSK. 

In  summary,  we  have  the  following  major  results  for  the  simulation  of 
the  narrow  bandpass-limited  satellite  channel: 

A.  The  best  modulations  studied  are  MSK  and  SQPSK.  At  a bit  error 
probability  of  ~ 0.005  a performance  loss  of  less  than  1 dB  and  3 dB  can  be 
expected  for  data  rates  up  to  32  kbps  and  48  kbps,  respectively,  through  a 
25-kHz  3-dB  bandwidth  FLEETSAT  channel,  if  either  modulation  scheme  is  em- 
ployed with  a coherent  matched-filter  receiver. 
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B.  The  next  best  modulations  are  QPSK  and  AQPSK.  QPSK  often  per- 
forms as  well  as  MSK  or  SQPSK,  and  any  additional  SNR  loss  for  QPSK  is  gen- 
erally less  than  1 dB.  AQPSK  performs  significantly  worse  than  MSK,  SQPSK 
and  QPSK  at  the  larger  data  rates,  but  can  still  yield  tolerable  losses  at 
reduced  rates. 

C.  CPQFSK  and  BPSK  are  the  worst  modulations  studied.  Although 
CPQFSK  does  much  better  than  BPSK,  neither  modulation  scheme  is  competitive 
with  the  others  at  the  larger  data  rates. 

D.  BPSK  would  suffer  a performance  loss  of  no  more  than  1 dB  at  a 
data  rate  of  19.2  kbps  for  a bit  error  probability  of  ~ 0.005  and  a 25-kHz 
3-dB  bandwidth  FLEETSAT  channel,  but  a loss  of  about  4 dB  at  a 24  kbps  rate. 

E.  For  a reliability  level  of  ^ 0.005  bit  error  probability  and  a 
performance  loss  of  1 dB  or  less,  the  maximum  standard  data  rate  possible 
through  a 25-kHz  3-dB  bandwidth  FLEETSAT  channel  is  19.2  kbps  for  BPSK, 

24  kbps  for  AQPSK  and  CPQFSK,  and  32  kbps  for  MSK,  SQPSK  and  QPSK. 

F.  BPSK  is  less  sensitive  to  an  offset  in  the  center  frequency  of 
the  FLEETSAT  channel  than  MSK  or  SQPSK.  (See  Fig.  3.4.) 

These  conclusions  are  generally  either  consistent  with  or  corroborated  by 
other  simulated  or  experimental  results  obtained  independently  [3-12].  But  a 
direct  comparison  of  our  data  with  that  of  the  literature  would  be  tenuous 
because  of  significant  differences  in  the  various  approaches.  Furthermore, 

such  a comparison  would  be  too  lengthy  for  inclusion  in  this  report. 

-2 

Error  rates  as  large  as  10  are  of  interest  in  coded  systems  whereas 

-3 

bit  error  probabilities  of  10  and  smaller  are  of  interest  in  uncoded  systems. 
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With  narrow  bandpass-limiting  operation  we  expect  larger  degradations  at 

the  lower  error  rates  since  there  the  limitations  are  more  likely  to  be 

filter  selectivity  and  intermodulation  effects  rather  than  channel  noise. 

-3 

Very  little  data  were  gathered  below  10  because  of  the  length  of  the 
simulations  required  with  both  uplink  and  downlink  noise.  Techniques  exist 
for  reducing  the  simulation  time  when  the  uplink  noise  is  negligible  [3,6]. 
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II.  S IMUL AT ION  MODEL 


The  bandlimited  channel  and  modulations  used  in  the  computer  simulation 
model  are  discussed  in  this  section. 

2.1  Channel  Model 

A block  diagram  of  the  simulation  model  was  shown  in  Fig.  1.1.  The  UHF 
satellite  model  consists  of  a narrow  bandpass  filter  of  nominal  3-dB  band- 
width B = 25  kHz  followed  by  an  ideal  bandpass  limiter  which  removes  the 
envelope  modulation  of  its  input.  Additive  white  Gaussian  noise  (AWGN)  is 
assumed  for  both  the  uplink  to  the  satellite  and  the  downlink  from  the  satel- 
lite. The  entire  simulation  is  performed  digitally  and  at  baseband  with 
complex  signals  and  noise.  The  coherent  performance  of  techniques  employing 
various  forms  of  antipodal  signaling  and  detection  filters  matched  to  the 
transmitted  signals,  and  operating  at  the  same  data  rate  R will  be  compared 
on  a common  basis. 

Data  rates  in  excess  of  the  filter  bandwidth  (B/R  < 1)  will  be  of  princi- 
pal interest.  Thus,  the  relatively  narrowband  filter  (generally  with  a non- 
linear phase  characteristic)  and  the  highly  nonlinear  limiter  will  tend  to 
constrain  the  signal  spectrum  and  distort  the  received  waveform  thereby 
introducing  a significant  amount  of  intersymbol  interference  (ISI) . Such  ISI 
will  be  the  principal  source  of  degradation  from  the  theoretical  matched- 
filter  performance  in  the  presence  of  AWGN.  For  moderate  uplink  and  downlink 

energy  contrast  ratios  of  interest  (E,/N  ) and  (E./N  ).  , respectively, 

b o up  b o down  J 

the  limiter  generally  has  a mitigating  effect  on  the  ISI.  That  is,  compared 
to  a linear  transponder  with  the  same  average  power,  the  hard-limiting  opera- 
tion can  actually  improve  overall  performance  in  some  cases. 
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2.1.1  Matched-Filter  Receiver 


The  uplink  E^/N^  is  defined  in  the  usual  way  as  the  ratio  of  the  useful 

signal  energy  per  data  bit  received  at  the  satellite  to  the  single-sided 
uplink  noise  power  spectral  density  [13].  The  downlink  E^/N^  is  defined 
similarly  assuming  that  all  the  satellite  transmitter  power  is  applied  to 
a useful  downlink  signal.  Of  course,  some  of  this  downlink  power  is  devoted 
to  ISI  produced  by  the  narrow  bandpass-limited  channel  and  to  repeating 
uplink  noise.  Hence,  not  all  the  downlink  signal  energy  is  really  useful 
for  detection,  especially  with  a receiver  matched  to  the  undistorted  uplink 
signals. 

It  has  been  recently  demonstrated  by  computer  simulation  [14]  that  a 1 
or  2 dB  improvement  in  performance  is  possible  with  a more  complex  receiver 
structure.  In  this  approach  the  matched  filters  are  replaced  by  a set  of 
crosscorrelators  operating  over  an  interval  spanning  most  of  the  ISI.  The 
reference  signals  are  generated  by  passing  a finite  set  of  transmitted  wave- 
forms through  an  equivalent  of  the  known  channel  in  the  absence  of  noise. 

This  procedure  approximates  maximum- likelihood  detection  but  is  feasible  only 
when  the  channel  is  known  and  time-invariant,  and  when  most  of  the  ISI  is 
confined  to  a few  signaling  intervals. 

The  maximum- likelihood  approach  using  the  Viterbi  algorithm  is  also  an 
interesting  possibility  [ 15 ] . A preliminary  investigation  for  BPSK  indicates 
that  nearly  all  the  residual  ISI  could  be  removed  with  a digital  processor 
of  modest  complexity  [15a]. 


\ 

\ 

t 


14 


Thus,  under  certain  conditions  it  is  feasible  to  improve  on  the 
(mis)  matched-filter  performance.  But  for  the  purpose  of  this  report  we 
simply  recognize  this  possibility  and  simulate  only  the  matched-filter 
model.  More  sophisticated  receivers  will  be  left  for  future  study. 

2.1.2  Ideal  Bandpass  Limiter 

Since  the  bandpass  limiter  can  be  implemented  with  a solid-state 
device  at  UHF  with  much  less  severe  AM  to  PM  conversion  problems  that 
can  plague  TWT  implementations  at  GHz  frequencies  [4,  5,  16],  the  ideal 
(hard-limiting)  bandpass-limiting  model  is  reasonable  for  our  simulation. 

The  ideal  bandpass-limiting  operation  is  depicted  in  Fig.  2.1  for  a 
baseband  complex  input  sample  x(n)  + jy(n),  marked  by  a clear  dot  in 
quadrant  I.  The  limiter  maps  this  point  to  the  circle  of  radius  r along 
a radial  from  the  origin  to  yield  a complex  output  sample  x* (n)  + jyT(n), 
marked  by  a solid  dot,  where 


Hence,  the  limiter  removes  the  amplitude  modulation  but  preserves  the 
phase  information,  i.e.,  for  a narrowband,  analog,  bandpass  input 


(2.1a) 


(2.1b) 


v(t)cos(oo  t + 4>(t))  at  a carrier  frequency  u)  / 2tt , the  ideal  bandpass- 
c c 
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IMAGINARY 


REAL 


Fig.  2.1.  Ideal  bandpass-limiting  operation. 


limiter  output  is  cosCoo^t  + <£(t)).  From  (2.1)  we  note  that  a change  in  either 
input  generally  produces  crosstalk  between  the  quadrature  outputs. 

By  assuming  that  the  limiter  leaves  the  signal-to-noise  ratio  (SNR) 
unchanged,  although  in  reality  the  SNR  is  altered,  and  ignoring  any  ISI 
effects,  the  approximate  effective  E^/N^  as  derived  in  Appendix  A is 


<vv 


ef  f 


(E,  /N  ) (E,  /N  ) , 

b o up  b o down 

(E.  /N  ) + (E,  /N  ),  + f 

d o up  b o down  R 


(2.2) 


where  (E^/Nq)  and  (E^/Nq)^ow^  are  the  uplink  and  downlink  SNR's  multiplied 
by  B/R,  respectively  [1,17].  If  the  SNR  at  the  input  is  large,  the  limiter 
doubles  the  SNR  [18,  p.  311,  Prob.  13].  However,  this  gain  in  signal  to  noise 
actually  cannot  be  exploited;  the  effective  end-to-end  ratio  might  as  well 
be  calculated  assuming  that  this  gain  is  nonexistent.  It  is  demonstrated 
in  Appendix  A that  this  (E^/N  ) ^ is  within  a decibel  of  the  true  E^/Nq,  ignor- 
ing any  ISI  effects  of  the  filter.  The  terminology  "balanced"  uplink  and  down- 
link simply  means  that  (E,  /N  ) = (E,  /N  ), 

b o up  bo  down. 

In  a system  using  FLEETSAT  the  links  will  usually  be  unbalanced  because 

of  different  terminal  capabilities  [1].  Note  from  (2.2)  that  if  one  of  the 

links  is  much  noisier  than  the  other,  then  (E,  /N  ) --  is  approximately  equal 

b o et t 

to  the  E^/N^  t*ie  noisY  link.  If  the  links  are  balanced,  then  (E^/NQ)e^^ 
is  comparable  to  and  at  most  3 dB  smaller  than  the  link  E^/N^’s.  In  this 
report,  B/R  will  usually  be  small  compared  to  the  sum  of  the  link  E^/N^’s. 
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Estimates  of  the  bit  error  probability  P^  obtained  by  counting  errors 


will  be  plotted  as  a function  of  (K  /N  ) --  in  Section  III.  For  the  sake  of 


comparison  solid  curves  representing  the  theoretical  bit  error  probability 


for  coherent  demodulation  of  equally-likely  antipodal  signals  in  AWGN  will 
also  be  shown.  The  degradation  from  this  theoretical  curve  due  to  the  narrow 
bandpass-limiting  channel  will  be  the  principal  quantity  of  interest. 

The  nonlinearity  introduced  by  the  bandpass  limiter  severely  hampers 
analysis  in  general.  Although  much  can  be  done  to  mitigate  ISI  in  a linear 
system  [15,  19],  and  though  bandpass  limiters  (without  filter  distortion) 
have  been  analyzed  to  a considerable  degree  [17,  20,  21],  the  narrow  band- 
pass-limited  channel  is  better  handled  through  simulation  or  experimentation. 
The  approaches  taken  in  much  of  the  literature  dealing  with  linear  channels 
even  resort  to  simulation  [3],  experimentation  [9],  or  restricted  analyses 
and  numerical  methods  [7,  8,  10,  11].  Since  the  satellite  components  were  un- 
available, we  chose  simulation. 

On  the  other  hand,  a promising  theoretical  approach  has  just  appeared 
[22]  for  handling  a bandlimited  nonlinear  channel  for  the  special  case  of 
BPSK  and  downlink  noise  only  using  the  Forney  scheme  [15].  Since  [22]  is  in 
the  realm  of  more  complex  receiver  structures  mentioned  earlier,  an  analysis 
of  this  sort  will  be  deferred.  Unfortunately,  it  will  probably  be  difficult 
to  extend  [22]  to  more  attractive  modulations  such  as  QPSK,  let  alone  MSK  or 
SQPSK. 


00 


(2.3) 
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We  know  of  just  three  other  investigations  [4-6]  that  are  closely  related 
to  our  objective  of  communi eating  at  high  data  rates  with  the  model  of  Fig. 
1.1.  However,  none  of  these  efforts,  which  were  experimental  in  nature, 
permitted  both  uplink  and  downlink  noise  sources  simultaneously  or  included 
all  the  modulations  we  consider. 

2.1.3  Digital  Processing 

Several  different  digital  filters  were  employed  in  the  simulation.  Fre- 
quency sampling  techniques  were  used  to  design  some  of  the  filters.  The  bi- 
linear transformation  from  the  analog  to  digital  domain  was  used  in  other 
design  procedures  for  the  Chebyshev,  Buttervorth  and  elliptic  filters  [23,24]. 
This  transformation  avoids  the  problem  of  aliasing  of  the  analog  frequency 
spectrum  in  the  digital  representation  of  a system  of  discrete  uniformly- 
spaced  samples.  But  the  effect  of  the  aliasing  on  the  accuracy  of  the  model 
is  small  if  the  analysis  bandwidth  W is  sufficiently  large. 

A baseband  simulation  involving  complex  signals  and  noise  is  more  con- 
venient mathematically  and  is  more  efficient  than  a bandpass  one  because  a 
much  smaller  sampling  rate  is  required.  The  channel  filter  3-dB  bandwidth  B 
should  be  distinguished  from  the  analysis  bandwidth  W,  the  complex  sampling 
rate  of  our  baseband  digital  simulation  model.  In  order  to  insure  a reason- 
ably accurate  approximation  to  the  bandpass  situation,  W must  be  signifi- 
cantly larger  than  the  greater  of  B and  the  3-dB  bandwidth  of  the  signaling 
spectrum  B , since  the  ideal  bandpass  limiter  generates  extensive  intermod- 
ulation products.  We  have  decided  on  the  rule-of-thumb 
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(2.4) 


W ¥ 4 max  {B,  B } 
s 


as  an  adequate  analysis  bandwidth. 

Some  experimentation  was  performed  to  test  the  sufficiency  of  this  rule- 
of-thumb.  We  found  that  doubling  complex  sampling  rates  approximately  equal 
to  the  right-hand  side  of  (2.4)  did  not  change  the  measured  error  probabilities 
appreciably.  Other  investigations  of  binary  matched-filter  detection  have 
concluded  that  four  complex  samples  per  data  bit  are  adequate  [25,  26].  Al- 
though this  conclusion  is  not  necessarily  valid  with  the  limiter  present  or 
when  Bg  is  significantly  larger  than  B,  it  is  consistent  with  (2.4),  since 
the  3-dB  bandwidth  of  BPSK  is  approximately  0.9R  (see  Section  2.2.2). 

2 . 2 Modulation  Schemes 

Only  modems  that  achieve  the  theoretical  matched-filter  performance  for 
antipodal  signals  and  the  distortionless  channel  in  AWGN  have  been  considered. 
The  modulations  to  be  simulated  are  denoted  by  the  following  acronyms  [27]: 


BPSK 


Binary  Phase  Shift  Keying 


QPSK 


QuadriPhase  Shift  Keying 


SQPSK 


Staggered  or  offset  QPSK;  same  as  QPSK  except 
that  phase  changes  are  restricted  to  0 or  ±tt/2 
and  may  occur  twice  as  often 


AQPSK 


Alternating  QPSK;  same  as  SQPSK  except  that  phase 
changes  are  restricted  to  ±tt/2 


20 


MSK 


Minimum  Shift  Keying;  continuous-phase  binary 
frequency  shift  keying  with  a frequency  separation 
of  R/2;  same  as  SQPSK  except  that  the  quadrature 
modulating  waveforms  are  sinusoidal  rather  than 
rectangular  in  shape 

CPQFSK  Continuous-Phase  QuadriFrequency  Shift  Keying  with 

a frequency  separation  of  R/4. 

Typical  modulation  waveforms  for  every  scheme  but  BPSK  are  drawn  in  Fig. 

2.2.  BPSK  is  like  one  channel  of  QPSK  but  with  half  the  intersymbol  spacing 
for  the  same  data  rate. 

2.2.1  Special  Properties 

Alternating  QPSK  is  of  special  interest  for  symbol  timing  acquisition 
and  tracking  because  a phase  transition  is  guaranteed  for  every  data  bit. 

AQPSK  also  has  a built-in  immunity  to  first-order  ISI  effects  [28]  . 

Staggered  QPSK,  MSK  and  AQPSK  are  distinguished  by  the  kindly  way  an 
ideal  bandpass  limiter  treats  their  filtered  waveforms.  Rhodes’  computer 
simulations  [29]  showed  that  SQPSK  and  MSK  retain  their  bandlimited  spectra 
with  low  sidelobes  after  ideal  bandpass  limiting,  in  contrast  to  QPSK.  Other 
experimentation  has  shown  that  the  limiter  builds  up  the  sidelobes  of  filtered 
AQPSK  but  to  magnitudes  well  below  their  prefiltered  levels  [28]  . 

When  there  is  little  or  no  sidelobe  build-up,  the  output  of  the  satellite 
power  amplifier,  which  the  bandpass  limiter  represents,  need  not  be  filtered 
as  much,  if  at  all.  Even  without  such  filtering,  however,  the  sidelobe  build- 
up with  QPSK  is  not  very  significant  in  terms  of  measured  error  probabilities 
as  will  be  seen  in  Section  III. 


21 


♦ 1 

COS  Ojt 

i 

CHANNEL  0 

2T 

4T 

sin  <*jt 

t 

CHANNEL  0 

2T 

4T 

a.  QPSK. 

cos  wt 

CHANNEL  0 

/ " 

ls 

v 

s. 

_ 

AT 

' msk/ 

\ 

OFFSET  QPSK 

sin  wt 

V' 

/ 

’N 

\ 

s. 

CHANNEL  0 

-( 

— 

T 

3T 

E 

\ 

~y 

5T 

b.  MSK  and  SQPSK. 


♦ 1 


COS  <jt L 

CHANNEL  0 


2T 


3T 


4 T 5T 


t 


♦ 1 

sin  ut 

CHANNEL  0 

• 1 

T 2T  3T 

4T  5T 

c.  AQPSK. 


+4 

cos  wt 

\ 1 

CHANNEL 

-1 

+1 

ZT 

* 

*in  wt 

1 \ 

J ^ 

CHANNEL  □ 

27  \ 

/at 

6T 

H 

d.  CPQFSK . 


Fig.  2.2.  Modulation  waveforms. 
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2.2.2  Signal  Bandwidths 

We  should  be  cognizant  of  the  relative  signaling  bandwidths  of  different 
modulations  when  interpreting  the  data  of  Section  III.  Let  the  double-sided 
r-dB  bandwidth  of  a baseband  signal  be  defined  as  twice  the  frequency  where 
the  power  spectral  density  first  falls  to  r dB  below  the  f=0  value.  And  let 
the  double-sided  zero-crossing  bandwidth  be  twice  the  frequency  of  the  first 
spectral  null.  Then  for  an  infinite  sequence  of  uniformly- distributed  symbols, 
several  of  these  bandwidths  are  listed  in  Table  2.1  for  each  modulation 
scheme. 

k 

We  are  assuming  without  verification  that  AQPSK  has  the  same  spectrum 

as  BPSK  [28],  and  that  SQPSK  has  the  same  spectrum  as  QPSK  [5,29],  Both  BPSK 

2 

and  QPSK  have  the  classical  (sin  r/r)  spectrum  with  a zero-crossing  bandwidth 
of  exactly  twice  the  symbol  rate  [30]. 

The  zero-crossing  bandwidth  for  MSK  is  exactly  3R/2,  and  the  MSK  3-dB 
bandwidth  of  Table  2.1  is  accurate  [3l]  . The  other  bandwidths  for  MSK  and 
all  the  CPQFSK  bandwidths  are  only  approximations.  The  CPQFSK  bandwidths  were 
taken  from  curves  in  [32]  . General  expressions  for  continuous-phase  FSK 
spectra  are  given  in  [32,  33  ] . 

From  Table  2.1  we  see  that  the  QPSK  bandwidth  is  half  that  of  BPSK,  as  is 
well  known,  and  that  the  r-dB  bandwidth  of  MSK  is  less  than  1.5  times  that  of 
QPSK.  However,  MSK  generally  has  a much  smaller  fraction  of  out-of-band  energy 
than  QPSK  [5  ] . The  spectrum  of  CPQFSK  is  also  quite  compact  because  of  the 
minimum  signaling  frequency  separation  for  orthogonality  and  the  phase 
continuity  [31,  32]  . 

k 

This  statement  may  be  an  assumption  of  fact. 
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TABLE  2.1 


DOUBLE-SIDED  SIGNALING  BANDWIDTHS  FOR  INFINITE  UNIFORMLY- 
DISTRIBUTED  DATA  SEQUENCES 


Bandwidth  (R  is  the  data  rate  In  bits  per  sec) 


Modulation 

3-dB 

13-dB 

20-dB 

Zero-Crossing 

BPSK,  AQPSK 

0.886R 

1.70  R 

1.92R 

2R 

QPSK,  SQPSK 

0.443R 

0.852R 

0.96R 

R 

MSK 

0.595R 

1.13  R 

1.32R 

3R/2 

CPQFSK 

0.63  R 

1.59  R 

1.62R 

^2.5R 
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III.  PERFORMANCE  ANALYSIS 


The  simulation  results  using  an  IBM  370/168  computer  are  presented  in 
this  section.  Unless  stated  otherwise  the  data  and  Gaussian  noise  sequences 
were  selected  using  the  nearly  uniformly  distributed  pseudorandom  number  gen- 
erator [34] 


where  is  a 9-digit  integer,  and  a polar  method  for  normal  deviates  [35]. 

* 

The  same  sequences  were  employed  when  comparing  different  modulations. 

Unless  stated  otherwise,  each  datum  shown  graphically  in  this  section 

represents  about  2000  bit  decisions;  some  runs  in  Fig.  3.4  were  made  with 

8000  trials  for  better  confidence  intervals.  Most  of  the  performance  estimates 

— 3 r^j  'v  — 2 

are  made  in  the  bit  error  probability  range  10  < P,  < 10  . Then  according 

b 

to  (2.3),  (E^/N^  ^^s  in  the  4 to  10  dB  range  are  of  primary  interest,  allow- 
ing for  a few  dB  degradation  from  the  theoretical  performance  curve. 

Typical  measurement  errors  for  a 95%  confidence  level  are  listed  in 
Table  3.1  for  2000  and  8000  bit  decisions  per  datum.  For  a given  measurement 


and  where  p specifies  the  confidence  interval.  Then  the  actual  loss  L is 

A A 

in  the  range  L < L < L+  with  probability  0.95. 

Some  digital  filters  used  in  the  simulations  are  listed  in  Table  3.2. 


Every  modulation  scheme  was  tested  in  the  presence  of  AWGN  only,  as  a check 
on  the  simulation  program;  performance  was  always  close  to  the  theoretical  bit 
error  probability  of  (2.3). 


K , - = ll5  K (mod  231  - 1) 
n+1  n 


(3.1) 


let  L and  L be  the  degradation  from  the  theoretical  curve  at 


+ 


(1  - p)P,  and  (1  + p)P_,  respectively,  where  P,  is  the  measured  error  rate 
b-  b b 


25 


TABLE  3.1 


RANGE  OF  ACTUAL  BIT  ERROR  PROBABILITY  FOR  95%  CONFIDENCE  LEVEL 


NUMBER  OF  BIT  DECISIONS 

A 

2000 

8000 

P. 

b 

P 

range 

P 

range 

O.OIO 

40% 

0.006  P.  ~ 0.014 

b 

20% 

0.008  < Pu  *2  0.012 
b 

0.005 

55% 

0.0022  ~ Pv  ~ 0.0078 
b 

29% 

0.0035  ~ P,  ~ 0.0065 
b 

0.001 

130% 

0 ~ P,  ~ 0.0023 

70% 

0.0003  < 0.0017 

* _ 

b 
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TABLE  3.2 


DIGITAL  FILTERS  SIMULATED 


Index 

Type 

Phase 

i 

Chebyshev 

Nonlinear 

2 

Butterworth 

Nonlinear 

3 

t 

FIR 

Linear 

4 

elliptic 

Nonlinear 

5 

FIR 

Nonlinear 

6 

FIR 

Cubic 

t Finite  Impulse  Response 
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The  Chebyshev,  Butterworth  and  elliptic  filters  were  implemented  in  the  time 
domain.  The  Finite  Impulse  Response  (FIR)  filters  were  realized  by  a fast 
convolution  method  via  the  frequency  domain.  Additional  information  is  given 
in  Appendix  B. 

For  programming  convenience  an  integral  number  of  complex  samples  per 
symbol  was  guaranteed  in  the  simulations.  An  adequate  analysis  bandwidth  was 
used  for  all  the  results  reported  in  this  section.  (See  Section  2.1.3.) 

The  principal  parameter  of  interest  for  a set  of  data  will  be  the  3-dB 
filter  bandwidth  to  data  rate  ratio  B/R.  Different  modulations  will  always 
be  compared  using  the  same  filter.  We  will  emphasize  results  for  the  filters 
that  come  closest  to  the  FLEETSAT  filter.  There  is  an  overall  performance 
comparison  among  different  filters  in  Section  I. 

3. 1 Linear-Phase  FIR  Filter 

Filter  3 of  Table  3.2  was  the  first  to  be  designed  [36]  and  simulated. 
This  linear  phase  filter  is  used  as  a baseline  for  comparison,  especially 
with  respect  to  results  for  nonlinear  phase  filters.  Its  frequency  response 
and  the  associated  performance  data  for  B/R  = 0.69  and  0.46  are  given  in  Figs. 
10-12  of  [37].  There  it  was  concluded  that: 

A.  A data  rate  of  32  kbps  is  feasible  with  little  or  no 
performance  degradation  for  QPSK,  SQPSK,  AQPSK,  or  MSK; 

B.  A data  rate  of  48  kbps  is  feasible  with  only  a slight 
increase  (<1  dB)  in  degradation  for  SQPSK  and  MSK;  the 


28 


degradation  increase  for  QPSK  and  AQPSK  is  approximately 
1 dB  and  2 dB,  respectively,  at  this  rate; 

C.  Neither  BPSK  nor  CPQFSK  is  feasible  with  acceptable 
performance  degradations  at  data  rates  of  32  kbps  or  above, 
although  CPQFSK  performed  significantly  better  than  BPSK; 

D.  For  bit-error  probabilities  down  to  10  ~ and  data  rates 
of  32  kbps  or  above,  the  best  performance  is  achieved  by 
MSK  followed  by  SQPSK,  QPSK,  AQPSK,  CPQFSK,  and  BrSK  in 
order  of  increasing  degradation. 

These  degradations  are  optimistic,  mainly  due  to  the  linear  phase  of 
filter  3,  as  we  shall  see  when  exploring  nonlinear-ohase  filters. 

The  predominance  of  errors  for  CPQFSK  resulted  from  signaling  with  the 
two  frequencies  at  the  band  edge,  where  more  signal  distortion  is  to  be 
expected.  The  inferior  performance  of  BPSK  is  attributable  to  the  facts 
that  the  signaling  bandwidth  of  BPSK  is  twice  that  of  QFSK  (see  Table  2.1) 
and  that  phase  transitions  are  always  of  180  degret  s with  BPSK  modulation. 

Rhodes  [29]  explains  how  polarity  changes  in  the  bandlimited  signal  with 
these  180-degree  phase  shifts  cause  destructive  ISI  and  a sharp  decrease  in 
the  signal  envelope.  The  hard-limiting  rtion  of  the  bandpass  limiter  tends 
to  restore  the  signal  envelope  to  its  prefiltered  value  while  maintaining  the 
rapid  phase  change.  Thus,  we  have  a spectral  sidelobe  build-up  with 
BPSK  as  with  QPSK  (see  Section  2.2.1). 

Although  AQPSK  apparently  has  the  same  bandwidth  as  BPSK,  there  is  less 
sidelobe  build-up,  evidently  because  AQPSK  has  only  90-degree  phase  changes. 
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This  may  help  explain  why  AQPSK  performs  better  than  BPSK  but  the  built-in 
immunity  of  AQPSK  to  ISI  mentioned  in  Section  2.2.1  is  a more  likely  explana- 
tion. BPSK  has  none  of  this  immunity. 

In  comparing  the  performance  of  SQPSK  and  QPSK,  which  apparently  have 
the  same  bandwidth,  we  see  from  conclusion  B above  that  the  sidelobe  build- 
up effect  is  relatively  minor  because  it  can  account  for  less  than  a 1-dB 
degradation  with  filter  3. 

Tests  between  fixed  and  pseudorandom  data  sequences  were  run  for  MSK, 
SQPSK,  AQPSK  and  QPSK  with  B/R  = 0.69.  For  fixed  sequence  "1"  the  modulating 
data  streams  alternated  between  +1  and  -1  on  both  quadrature  channels;  for 
sequence  M2"  one  data  stream  alternated  while  the  quadrature  stream  was  held 
constant.  These  sequences  were  likely  candidates  for  generating  a maximal 
amount  of  ISI  and/or  crosstalk  (see  Section  2.1.2).  The  crosstalk  phe- 
nomenon due  to  the  bandpass  limiter  has  been  verified  experimentally  [5,  38]. 

With  MSK  and  filter  3,  sequence  1 yielded  performance  a fraction  of  a 
decible  worse  than  a pseudorandom  data  stream,  but  sequence  2 did  slightly 
better  than  a pseudorandom  sequence.  With  SQPSK  sequence  2 performed  about 
1 dB  worse  than  sequence  1 but  not  any  worse  than  a pseudorandom  sequence. 
With  AQPSK  sequence  2 was  slightly  worse  than  1 but  neither  sequence  appeared 
to  be  as  bad  as  a pseudorandom  sequence.  Thus,  these  contrived  sequences  were 
not  harmful  to  the  offset  modulations.  However,  sequence  1 yielded  a 1-dB 
degradation  over  pseudorandom  performance  with  QPSK. 

3 . 2 Butterworth  and  Chebyshev  Filters 

With  nonlinear-phase  Butterworth  and  Chebyshev  filters,  we  found  that 
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sequences  1 and  2 were  far  from  the  worst-case  inputs  for  MSK  and  AQPSK, 
respectively.  As  shown  in  Fig.  3.1,  pseudorandom  sequences  gave  many  more 
errors  than  these  fixed  sequences.  Therefore,  it  is  not  obvious  what  the 
worst-case  data  sequence  would  be  for  a given  offset  modulation  and  satellite 
model. 

Note  from  Fig.  3.1  that  AQPSK  performance  is  approximately  2 dB  and 
2.5  dB  worse  than  MSK  for  the  Butterworth  and  Chebyshev  filters,  respectively. 
Also  observe  that  the  losses  for  the  Chebyshev  filter  are  about  1 dB  wo’\se 
than  the  Butterworth  filter,  even  though  the  data  rate  is  13%  smaller.  As 
explained  in  Section  I the  Butterworth  filter  characteristics  (see  Appendix  B) 
are  closer  to  those  of  the  FLEETSAT  filter. 

-2 

In  Fig.  3.1a  we  see  that  MSK  suffers  a 2-dB  degradation  for  P^  « 10 
and  B/R  = 0.54.  Comparing  this  with  the  filter-3  data  (see  Fig.  1.3b),  the 
nonlinear  phase  of  the  Butterworth  filter  is  probably  responsible  for  no  more 
than  1 dB  of  this  performance  loss,  with  the  slightly  faster  spectral  roll-off 
accounting  for  the  rest. 

Additional  data  gathered  under  the  conditions  of  Fig.  3.1  except  for 
unbalanced  links  was  not  sufficiently  different  from  the  balanced  case  to 
warrant  graphical  display.  For  a noisy  uplink  and  good  downlink  ((E^/Nq)^^ 

= 20  dB)  we  observed  about  1 dB  of  additional  degradation  compared  to  the 
balanced  case.  With  a good  uplink  ((E^/N^)  = 20  dB)  and  a noisy  downlink, 

performance  was  about  the  same  as  in  the  balanced  case  except  with  MSK  and  the 
Butterworth  filter,  where  an  additional  less  of  about  1 dB  was  observed. 


a.  Butterworth  filter;  B/R  = 0.54. 


Fig.  3.1.  MSK  and  AQPSK  performance  for  Butterworth  and  Chebyshev  filters, 
bandpass  limiter,  random  and  non-random  data,  and  balanced  links. 
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0.1 


b.  Chebyshev  filter;  B/R  = 0.62. 


Fig.  3.1.  Continued. 


3.3  Elliptic  Filter 


An  elliptic  filter  with  about  the  same  spectral  roll-off  but  a much  more 
nonlinear  phase  characteristic  (see  Appendix  B)  was  simulated  next.  By  com- 
paring the  MSK  and  AQPSK  performance  of  Fig.  3.2a  with  that  of  Fig.  3.1, 
we  see  that  this  filter  represents  a pessimistic  extreme  in  the  same  way  that 
the  linear-phase  FIR  filter  represented  an  optimistic  extreme.  As  explained 
in  Section  I the  Butterworth  filter  and  the  12-degree  version  of  filter  6 to 
be  discussed  later  are  closer  to  the  FLEETSAT  filter.  Nevertheless,  we  include 
the  data  of  Fig.  3.2  as  an  interesting  comparison  to  the  Butterworth  and 
Chebyshev  simulations. 

At  3/4  the  data  rate  the  performance  loss  for  MSK  in  Fig.  3.2a  is  reduced 
to  2 dB  in  Fig.  3.2b  but  the  loss  for  AQPSK  is  still  unacceptable.  This  is  the 
only  filter  simulated  for  which  the  MSK  performance  was  significantly  worse 
than  any  of  the  other  modulations. 

The  data  of  Fig.  3.2a  for  MSK  and  AQPSK  were  basically  unchanged  by  a 
four-fold  increase  in  the  sampling  rate  using  an  elliptic  filter  with  other- 
wise identical  characteristics. 

With  fixed  data  sequence  1 and  the  conditions  of  Fig.  3.2a,  QPSK  yielded 
just  a few  more  errors  than  a pseudorandom  sequence.  Sequences  1 and  2 both 
resulted  in  less  than  a 1-dB  loss  for  MSK,  as  did  sequence  1 for  SQPSK  and 
AQPSK;  for  SQPSK  sequence  2 performed  worse  but  still  significantly  better 
than  a pseudorandom  stream,  and  with  only  a 1.5-dB  loss  for  AQPSK.  For  the 
smaller  data  rate  of  Fig.  3.2b,  sequence  1 for  MSK  and  sequence  2 for  SQPSK 
and  AQPSK  performed  better  than  but  much  closer  to  a pseudorandom  sequence, 
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a.  B/R  = 0.72. 

Fig.  3.2.  Performance  for  elliptic  filter,  bandpass  limiter,  random  data,  and 
balanced  links. 
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b.  B/R  = 0.96. 


Fig.  3.  2.  Continued. 


36 


\ 


as  would  be  expected.  Thus,  neither  fixed  sequence  is  a worst-case  data 
stream  for  offset  modulations  and  our  satellite  model.  Hence,  all  simulations 
in  the  sequel  employed  pseudorandom  data. 

3. 4 BPSK  Performance 

We  must  drastically  reduce  the  data  rate  with  BPSK  to  get  acceptable 
performance  losses  as  shown  in  Fig.  3.3.  With  the  exception  of  the  filter-3 
data  which  shifted  to  the  right  about  1 dB,  the  performance  for  a solid  up- 
link ((E^/NQ)up  = 20  dB)  is  close  to  that  of  Fig.  3.3a.  In  these  cases  the 
faster  spectral  roll-off  and  nonlinear  phase  of  the  Butterworth  filter  are 
each  responsible  for  about  a 1-dB  degradation  over  filter  3,  taking  discrep- 
ancies in  B/R  ratios  into  account.  This  is  consistent  with  earlier  MSK 
results. 

In  Fig.  3.3b,  where  the  uplink  E^/N^  varies  between  5 and  7 dB,  the 
Chebyshev  filter  degrades  almost  1 dB  more  than  in  the  balanced  case.  For 
the  noisy  uplink  case  of  Fig.  3.3c,  the  performance  for  filter  3 and  the 
Butterworth  filter  is  about  1 dB  worse  than  in  Figs.  3.3a  and  b;  the 
Chebyshev  and  elliptic  filters  do  only  a fraction  of  a decibel  worse  than 
before . 

We  believe  that  the  noticeably  smaller  losses  with  balanced  links  com- 
pared to  a noisy  uplink  case  are  because  the  bandpass  limiter  mitigates  the 
ISI  introduced  by  the  filter  when  the  uplink  is  not  dominated  by  noise. 

3 . 5 Nonlinear-Phase  FIR  Filter 

The  next  model  simulated,  filter  5 of  Table  3.2,  has  characteristics  very 
close  to  the  FLEETSAT  filter  (see  Fig.  1.2).  In  an  earlier  version  the  phase 
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a.  Balanced  links. 


Fig.  3.3.  BPSK  performance  for  various  filters,  bandpass  limiter,  and  random 
data. 
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Fig.  3.3.  Continued. 
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Fig.  3.3.  Continued. 
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response  at  the  band  edges  was  specified  rather  arbitrarily  in  the  absence  of 
available  phase  deviation  data  and  is  probably  too  nonlinear.  In  this  respect 
the  filter-5  simulation  data  presented  in  [39]  was  conservative. 

The  performance  of  the  latest  FLEETSAT  filter  model  (filter  5)  is  shown 
in  Fig.  3.4  for  BPSK,  MSK  and  SQPSK.  Most  of  the  points  represent  8000  bit 
decisions.  Figures  3.4a  and  b present  data  for  BPSK  and  MSK  and  balanced  links 
where  (E^/N^)  = ^^b^o^down’  ^.4c  gives  data  for  unbalanced  links.  The 

performance  of  SQPSK  was  slightly  (but  not  significantly)  worse  than  that  of 
MSK  in  Fig.  lc  and  with  balanced  links,  based  on  data  not  shown. 

For  a zero  center  frequency  offset  there  is  less  than  a one  decibel  degra- 
dation in  BPSK  for  B/R  = 1.35.  However,  the  loss  increases  rapidly  with  the 
data  rate  to  about  4 dB  at  P^  = 0.005  for  B/R  = 0.94.  With  MSK  there  is  very 
little  loss  for  B/R  = 0.81  but  a 2 or  3-dB  degradation  for  B/R  = 0.54. 

Additional  performance  losses  due  to  an  uncompensated  center  frequency 
offset  in  the  FLEETSAT  filter  are  also  indicated  in  Figs.  3.4a  and  b.  A 2% 
frequency  offset  relative  to  the  3-dB  filter  1 andwidth  B at  B/R  = 1.35  seems 
tolerable  for  BPSK.  But  higher  percentage  offsets  mean  much  larger  losses 
as  shown  in  Fig.  3.4a.  Even  a 1%  frequency  offset  at  B/R  « 0.81  and  0.54  means 
a significant  additional  loss  for  MSK  and  SQPSK,  however.  Keeping  in  mind  the 
signaling  bandwidth  differences  for  these  modulations  (see  Table  2.1),  we  con- 
clude that  MSK  and  SQPSK  are  more  sensitive  than  BPSK  to  uncompensated  center 
frequency  offsets  in  the  FLEETSAT  channel. 

Figure  3.4c  demonstrates  that  the  SNR  enhancement  of  the  limiter  for  a 
solid  uplink  cannot  be  exploited,  since  the  performance  loss  is  no  less  than 
with  balanced  links  or  a solid  downlink.  (S:*e  Section  2.1.2  and  Appendix  A.) 


41 


a.  BPSK;  balanced  links. 


Fig.  3.4.  Performance  for  FIR  filter  5 (closest  filter  to  FLEETSAT) , bandpass 
limiter,  random  data  and  balanced  links. 
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b.  MSK;  balanced  links. 


Fig.  3.4.  Continued. 
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c.  MSK  and  SQPSK;  unbalanced  links. 


Fig.  3.4.  Continued. 
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3.6  Cubic-Phase  FIR  Filter 


The  last  filter  simulated,  FIR  filter  6 of  Table  3.2  and  Appendix  B, 
has  an  overall  magnitude  response  that  is  quite  faithful  to  that  expected  for 
the  FLEETSAT  filter  and  an  adjustable  linear,  quadratic  or  cubic-phase  char- 
acteristic. In  the  nonlinear-phase  cases,  an  arbitrary  phase  deviation  from 
linearity  can  be  chosen. 

Filter  6' s performance  for  B/R  = 0.5  and  all  modulations  but  BPSK  is 
shown  in  Fig.  3.5.  In  Figs.  3.5a  and  b we  compare  results  for  12  and  36- 
degree  phase  deviations  at  0.35B  to  the  right  of  band  center.  Figures  3.5b 
and  c compare  performance  for  cubic  and  quadratic  phases.  Finally,  Figs.  3.5b 
and  d show  results  for  different  correlation  delays. 

MSK  yields  the  best  performance  while  CPQFSK  does  very  badly.  In  con- 
trast to  the  filter-3  data,  there  was  no  predominance  of  errors  when  signaling 
at  the  two  extreme  frequencies. 

Performance  is  rather  insensitive  to  a change  in  correlation  delay  of 
25%  of  the  integration  time  of  2/R  for  MSK,  SQPSK,  QPSK  and  CPQFSK  [27]; 
previously,  a 12.5%  change  for  MSK  made  little  difference  with  the 
Chebyshev  filter  [38].  It  is  remarkable  that  a change  in  the  AQPSK  correlation 
delay  of  50%  of  the  integration  time  of  1/R  was  not  more  noticeable,  but  the 
greater  tolerance  for  first-order  ISI  would  be  a contributing  factor. 

For  reasons  discussed  in  Section  I we  predict  a 2 or  3-dB  performance 
loss  for  B/R  = 0.52  MSK  or  SQPSK  and  the  FLEETSAT  model  based  largely  on  the 
12-degree  filter-6,  Butterworth  and  filter-5  simulations. 

Worst-case  filter-6  performance  for  BPSK  at  various  B/R  ratios  is  shown 
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a.  Cubic  phase,  12°  deviation,  15.7/R  delay. 


Fig.  3.5.  Performance  for  FIR  filter  6,  various  modulations,  bandpass  limiter, 
random  data  and  balanced  links  (B/R  = 0.505,  deviation  at  0.35B). 
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b.  Cubic  phase,  36°  deviation,  15.7/R  delay. 


Fig.  3,5.  Continued. 
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c.  Quadratic  phase,  36°  deviation,  15.7/R  delay. 


Fig.  3.5.  Continued. 
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d.  Cubic  phase,  36°  deviation,  15.2/R  delay. 


Fig.  3.5.  Continued. 
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in  Fig.  3.6.  In  view  of  the  sharper  roll-off  of  the  FLEETSAT  filter  these 
losses  are  somewhat  optimistic  at  the  larger  data  rates  (see  Fig.  1.3a). 

Most  of  the  filter-6  simulation  data  can  be  condensed  into  Fig.  3.7 
which  gives  the  approximate  degradation  as  a function  of  the  data  rate 
to  filter  3-dB  bandwidth  ratio  R/B.  The  high  degradation  portions  of  the 
BPSK  and  CPQFSK  curves  were  based  on  loss  data  considerably  above  4 dB. 

Based  on  Fig.  3.7  the  maximum  feasible  data  rates  for  each  modulation 
scheme  are  listed  in  Table  3.3.  These  estimates  turned  out  to  be  somewhat 
optimistic  when  considering  the  greater  selectivity  of  the  FLEETSAT  filter. 
But  as  before,  we  observe  that  MSK  and  SQPSK  are  the  best  modulations  for  the 
narrow  bandlimited  channel,  followed  by  QPSK,  AQPSK,  CPQFSK  and  BPSK. 
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Fig.  3.6.  BPSK  performance  for  cubic-phase  filter,  bandpass  limiter,  random 
data,  and  balanced  links  (36°  deviation;  15.2/R  delay). 
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PERFORMANCE  LOSS  (dB) 


R/B 


a.  12°  deviation,  15.7/R  delay. 


Fig.  3.7.  Estimated  degradation  from  theoretical  at  P 0.005  vs.  R/B  for  cubic- 
phase  filter  6,  bandpass  limiter,  random  data,  and  balanced  links. 
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b.  36°  deviation,  15.2/R  delay. 


Fig.  3.7.  Continued. 
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TABLE  3.3 


* 

ESTIMATED  MAXIMUM  FEASIBLE  DATA  RATES 

THROUGH  THE  25-kHz  SATELLITE  MODEL  WITH  FILTER  6 

FOR  A 2-dB  PERFORMANCE  LOSS  AT  PL  « 0.005 

b 


Modulation 

Cubic  Phase  Deviation  at 
12°  8.8  kHz  from  Band  Center 

36° 

BPSK 

28  kbps 

20  kbps 

MSK 

56  kbps 

40  kbps 

SQPSK 

56  kbps 

40  kbps 

AQPSK 

48  kbps 

36  kbps 

QPSK 

52  kbps 

36  kbps 

CPQFSK 

44  kbps 

32  kbps 

given  to  the  nearest  A kbps 
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APPENDIX  A 


Derivation  of  (En  /N  ) 

d o eff 

Let  B be  the  bandwidth,  let  P and  P,  be  the  received  signal  power 

up  down 

on  the  uplink  and  downlink,  respectively,  and  let  (N  ) and  (N  ) , be  the 

o up  o down 

corresponding  single-sided  AWGN  power  spectral  densities.  Suppose  that  we 
can  ignore  any  ISI  effects  and  that  the  ideal  bandpass  limiter  alters  the  in- 
put SNR  by  a factor  y,  i.e.,  suppose  that  the  output  SNR  is 


SNR 

out 


ySNR 


in 


(A.l) 


where  we  write  (P/N  ) for  P / (N  ) . Then  the  fraction  of  received  down- 

o up  up  o up 

link  signal  power  devoted  to  the  uplink  signal  is  yP  / (yP  + (N  ) B)  and 

up  up  o up 

the  fraction  devoted  to  uplink  noise  is  (N  ) B/ (yP  + (N  ) B) . The  effective 

o up  up  o up 

SNR  at  the  receiver  can  be  defined  as  the  ratio  of  the  received  signal  power  de- 


voted to  uplink  signal  and  the  sum  of  the  downlink  noise  power  (N  ) , B and 

o down 

the  received  signal  power  devoted  to  uplink  noise,  i.e., 

yp. 


SNR 


(P/Veff 


( 


_yE_ 


eff 


yP  + (N  ) B 
up o up 


-)P 


down 


(N  ) 
o down 


(N  ) B 

B + ( SJffi 


(A. 2) 


r)P 


yP  + (N  ) B down 
up  o up 


Since  the  power  P can  be  generically  expressed  in  terms  of  the  energy 
per  bit  E^  and  the  data  rate  R as  P = E^R,  we  can  rewrite  (A. 2)  as 
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(A.  3) 


Y(E.  /N  ) (E  /N  ), 

(E./N  ) — — P~  Rb  ° dOWn • 

y(E  /N  ) + £ + (K  /N  ) 

b o up  R do  down 

Now  (2.2)  follows  by  setting  the  y of  (A.l)  equal  to  unity  in  (A. 3). 

We  now  argue  that  y=l  results  in  a good  approximation  to  the  true  E^/N^  by 
considering  the  limiting  cases  of  strong  and  weak  links,  where  we  take 
E^/N^  >>  1 ^°r  a stronS  and  E^/N^  <<:  1 f°r  a weak  link.  We  further 

assume  that  B/R  ~ 1. 

Suppose  first  that  .the  uplink  is  strong,  i.e.,  (E^/N^)^  >:>  !•  As  men“ 
tioned  in  Section  2.1.2  this  implies  that  y=2.  If  the  downlink  is  weak,  i.e., 

if  (Eb/No)down  <<  then  (A'3)  reduces  to  <Weff  ~ Wdown’  indePendent 
of  y,  which  agrees  exactly  with  our  y=l  approximation  (2.2)  in  this  downlink 

noise  limited  case.  If  the  downlink  is  strong,  there  are  two  subcases  to  con- 
sider. In  the  first  subcase,  (K  /N  ) » (E,  /N  ) , >>  1,  and  again 

d o up  b o down 

(E^/N^)^^^  « (E^/N^)^^^  from  both  (A. 3)  and  (2.2);  we  are  still  downlink  noise 
limited. 

In  the  second  subcase,  (E./N  ),  » (E  /N  ) >>  1,  and  (E  /N  ) -- ~ 

9 b o down  b o up  b o eff 

y(E,  /N  ) = 2(E,  /N  ) . In  this  instance  the  model  leading  to  (A. 3)  is  faulty 

d o up  d o up 

because  the  true  effective  E^/N^  cannot  exceed  the  (E^/N^)  of  the  weakest 
link,  i.e. , the  system  is  not  Gaussian.  This  has  been  verified  experimentally 
in  Fig.  3.4c,  for  example.  With  y=l,  our  approximation  (2.2)  is  okay  if  the 
ISI  effects  are  not  too  great,  or  provided  we  process  sufficiently  at  the  re- 
ceiver to  mitigate  the  ISI. 
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(A.  3) 


Y(E  /N  ) (E  /N  ), 

(E./N  ) b o up  b ° down . 

>(VVup  + R + <VNo>dovn 

Now  (2.2)  follows  by  setting  the  y of  (A.l)  equal  to  unity  in  (A. 3). 

We  now  argue  that  y=l  results  in  a good  approximation  to  the  true  E./N  by 

b o 

considering  the  limiting  cases  of  strong  and  weak  links,  where  we  take 
E^/N^  >:>  1 f°r  a stronS  li-nk  and  E^/N^  <<  1 for  a weak  link.  We  further 
assume  that  B/R  1. 

Suppose  first  that  the  uplink  is  strong,  i.e.,  (E^/N^)^  >>  !•  As 

tioned  in  Section  2.1.2  this  implies  that  y=2.  If  the  downlink  is  weak,  i.e., 

if  (E./N  ),  <<  1,  then  (A. 3)  reduces  to  (E,  /N  ) -_  «,  (E./N  ),  , independent 

b o down  b o eff  bo  down 

of  y,  which  agrees  exactly  with  our  y=l  approximation  (2.2)  in  this  downlink 
noise  limited  case.  If  the  downlink  is  strong,  there  are  two  subcases  to  con- 
sider. In  the  first  subcase,  (E,  /N  ) » (E,  /N  ) , » 1,  and  again 

b o up  bo  down 

(E^/NQ)e^^  ~ ^^b^o^down  trom  both  (A.  3)  and  (2.2);  we  are  still  downlink  noise 
limited. 

In  the  second  subcase,  (E  /N  ) , » (E  /N  ) » 1,  and  (E  /N  ) ^ ~ 

b o down  b o up  b o eff 

y(E^/NQ)u^  **  2(E^/No)u^.  In  this  instance  the  model  leading  to  (A. 3)  is  faulty 
because  the  true  effective  E^/N^  cannot  exceed  the  (E^/Nq)  of  the  weakest 
link,  the  system  is  not  Gaussian.  This  has  been  verified  experimentally 
in  Fig.  3.4c,  for  example.  With  y=l,  our  approximation  (2.2)  is  okay  if  the 
ISI  effects  are  not  too  great,  or  provided  we  process  sufficiently  at  the  re- 
ceiver to  mitigate  the  ISI. 
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Now  suppose  that  the  uplink  is  weak  so  that  y=7t/4  [18,  p.  311,  Prob.  13]. 


If  the  downlink  is  strong,  (A.  3)  reduces  to  (E^/N^)^ 

f ^Wup 

about  1 dB  smaller  than  the  (E,  /N  ) obtained  with  Y=1  from  (2.2) 

d o up 

if  the  downlink  is  weak,  (A. 3)  reduces  to  « (yR/B) (E^/Nq 

which  again  is  about  1 dB  smaller  than  the  y=l  approximation. 


which  is  only 
Finally, 

upWdown’ 
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APPENDIX  B 


Filter  Characteristics 

Some  properties  of  the  principal  filters  simulated  are  listed  in  Table 
B.l  using  standard  digital  signal  processing  terminology  [40].  We  will  pro- 
vide the  magnitude  and  phase  responses  for  filters  1,  2,  4 and  6 and  the 
impulse  responses  of  filter  6 in  this  appendix.  Frequency  responses  for 
filters  3 and  5 were  reported  in  [37]  and  [39],  respectively. 

The  Chebyshev  and  Butterworth  filters  were  designed  using  [41]  and 
realized  in  parallel  canonical  form,  whereas  the  elliptic  filter  was  obtained 
utilizing  a readily  available  design  program  [42]  and  implemented  in  cascade 
canonical  form.  Hence,  these  causal,  Infinite  Impulse  Response  (HR)  filters 
were  realized  entirely  in  the  time  domain. 

The  computed  responses  for  the  IIR  filters  are  shown  in  Figs.  B.l  and 
B.2.  The  Chebyshev  filter  is  equiripple  in  the  passband,  the  Butterworth  fil- 
ter is  maximally  flat  at  f=0,  and  the  elliptic  filter  is  equiripple  in  both  the 
passband  and  the  stopband.  These  magnitude  properties  are  preserved  from  the 
analog  domain,  thanks  to  a pre-warping  of  the  critical  analog  frequencies  and 
the  use  of  the  bilinear  transformation  during  the  design  process  [23,  24]. 

Using  frequency  sampling  techniques,  filters  3 and  6 were  selected  from 
[36]  and  filter  5 was  modeled  according  to  the  FLEETSAT  filter  characteris- 
tics of  Section  1.1.  All  three  FIR  filters  were  simulated  nonrecursively  by 
the  overlap-add  method  [23,  24].  This  involved  sectioning  the  input  sequence 
into  equal  segments,  generally  much  shorter  than  the  impulse  response  of  the 
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DIGITAL  FILTERS  SIMULATED 
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Fig.  B.l.  Response  of  Chebyshev  and  Butterworth  filters. 
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Fig.  B.2.  Response  of  elliptic  filter. 
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filter,  padding  with  complex  zeros  at  the  end  of  each  segment  to  twice  the 
response  length,  and  taking  fast  Fourier  transforms  (FFT’s).  Multiplying  the 
frequency  transforms  of  appropriately  padded  input  and  impulse  response  se- 
quences corresponds  to  a convolution  of  the  two  finite  length  sequences. 
Naturally,  the  impulse  response  transform  need  be  found  only  once. 

Although  every  FIR  filter  can  be  realized  with  a frequency  sampling 
structure  in  the  time  domain,  they  are  usually  implemented  with  the  above 
fast  convolution  technique.  However,  a time  domain  realization  can  require 
fewer  multiplications  if  nearly  all  the  frequency  samples  are  zero  [43]. 

The  impulse  responses  for  all  but  filter  3 were  non-anticipatory , 
i.e.,  causal.  Valid  results  can  be  obtained  by  simulating  non-causal  filters, 
however,  the  only  differences  being  in  the  output  delay  and  programming 
convenience . 

With  filter  3 we  padded  in  the  middle  of  the  impulse  response  since  the 
linear-phase  filter  designs  of  [36]  yield  impulse  responses  in  a non-causal 
form  with  the  peak  at  the  origin.  Alternatively,  we  could  have  rotated  the 
impulse  response  by  half  its  length  to  obtain  the  causal  form  and  then  padded 
with  an  equal  number  of  zeros  on  either  side  as  we  did  with  filter  6.  We 
padded  at  the  end  of  the  impulse  response  with  filter  5.  In  each  case  the 
number  of  complex  zeros  inserted  equaled  the  impulse  response  length,  which 
means  the  frequency  response  was  interpolated  to  twice  the  number  of  original 
frequency  samples. 
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A causal  linear-phase  FIR  filter  with  a real  impulse  response  has  a phase 
delay  of  (N-l)/2  samples  and  an  even  impulse  response,  i.e.,  h(N-l-n)  = h(n), 
where  h(n)  is  the  nth  sample  of  a response  of  length  N [24,  p.  78].  The 
main  worry  with  the  nonlinear-phase  FIR  filters  5 and  6 is  the  fact  that  their 
impulse  responses  can  be  rather  uneven,  as  we  shall  see  for  filter  6.  We  have 
verified  that  the  responses  for  filters  5 and  6 have  negligible  imaginary  parts, 
so  one  expects  the  peak  of  the  causal  impulse  response  to  appear  at  about  the 
(N/2)th  sample,  as  in  the  linear-phase  case. 

For  all  filters  simulated  the  correlation  delay  in  seconds  is  related 
to  the  number  of  samples  in  the  delay  by 

Td  = Nd/W  (A. 1) 

where  W is  the  complex  sampling  rate. 

The  impulse  and  magnitude  responses  for  the  cubic-phase  filter  6 are 

shown  in  Fig.  B.3.  The  phase  deviation  from  linearity  in  degrees  was  speci- 
4 3 

fied  by  D(f/10  ) , where  D is  either  12  or  36  degrees  and  where  f is  the  fre- 
quency deviation.  In  superimposing  this  phase  deviation  on  a causally  rotated 
linear-phase  design  selected  from  [36],  we  see  from  Fig.  B.3a  that  the  impulse 
response  has  become  asymmetrical  about  63/2  and  that  the  peak  has  shifted  to 
the  left  for  the  36-degree  case,  but  only  by  2 samples.  These  impulse  responses 
were  padded  on  each  side  with  32  zeros  in  obtaining  the  interpolated  frequency 
response.  The  correlation  delays  of  Table  B.l  for  filter  6 resulted  from  the 
sum  of  the  32  samples  on  the  left  and  the  location  of  the  peaks  in  Fig.  B.3a. 


64 


a.  Impulse  response. 


Fig.  B.3.  Characteristics  of  cubic-phase  FIR  filter  6. 
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Fig.  B.3.  Continued. 
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We  verified  that  changing  the  sign  of  the  cubic-phase  deviation  results 
in  an  impulse  response  that  is  essentially  a reflection  of  the  response  of 
Fig.  B. 3a  about  a vertical  line  through  its  peak.  In  other  words,  the  larger 
sidelobes  appear  after  the  main  lobe,  as  with  the  HR  filter  impulse  responses 
which  all  had  negative  phase  deviations  (see  Figs.  B.l  and  B.2).  By  this 
symmetry,  the  sign  of  the  phase  deviation  should  not  affect  the  simulated 
results  for  filter  6. 

A 2 

A quadratic-phase  version  with  the  phase  deviation  of  sgn(f)  . D(f/10  ) 
was  also  simulated  (see  Fig.  3.5c).  In  this  instance  the  maximum  sidelobe 
level  was  -78  and  -58  dB  for  D=12°  and  36°,  respectively  (see  Fig.  B.3b). 

The  asymmetry  in  the  impulse  responses  for  this  quadratic  case  was  less  pro- 
nounced than  in  Fig.  B.3a. 
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GLOSSARY 


Abbreviations 


AM 

Amplitude  Modulation 

AQPSK 

Alternating  QPSK 

AWGN 

Additive  White  Gaussian  Noise 

BPSK 

Binary  Phase  Shift  Keying 

CPBFSK 

Continuous-Phase  Binary  Frequency  Shift  Keying 

CPQFSK 

Continuous-Phase  QuadriFrequency  Shift  Keying 

FFT 

Fast  Fourier  Transform 

FIR 

Finite  Impulse  Response 

FLEETSAT 

unofficial  designation  for  a military  communications  satellite 

HR 

Infinite  Impulse  Response 

ISI 

InterSymbol  Interference 

MSK 

Minimum  Shift  Keying 

PM 

Phase  Modulation 

QPSK 

QuadriPhase  Shift  Keying 

SNR 

Signal- to-Noise  Ratio 

SQPSK 

Staggered  (offset)  QPSK 

TWT 

Traveling  Wave  Tube 

UHF 

Ultra  High  Frequency 

VHF 

Very  High  Frequency 

WSC-3 

satellite  communications  set  AN/WSC-3 
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Symbols 


B 

satellite  filter  3-dB  bandwidth 

received  signal  energy  per  data  bit 

<Eb/No>do»n 

downlink  E^/N^ 

<VVeff 

effective  E^/N^ 

<VV„p 

uplink  E^/N^ 

f 

frequency  deviation 

N 

o 

single-sided  noise  power  spectral  density 

pb 

bit  error  probability 

R 

data  rate 

W 

analysis  bandwidth  (complex  sampling  rate) 
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